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Percolation Thresholds of the Fortuin-Kasteleyn Cluster for a Potts 
Gauge Glass Model on Complex Networks 

Analytical Results on the Nishimori Line 

Chiaki Yamaguchi 
Kosugichou 1-359, Kawasaki 211-0063, Japan 

It was pointed out by de Arcangelis et al. [Europhys. Lett. 14 (1991), 515] that the 
correct understanding of the percolation phenomenon of the Fortuin-Kasteleyn cluster in the 
Edwards-Anderson model is important since a dynamical transition, which is characterized 
by a parameter called the Hamming distance or damage, and the percolation transition are 
related to a transition for a signal propagating between spins. We show analytically the 
percolation thresholds of the Fortuin-Kasteleyn cluster for a Potts gauge glass model, which 
is an extended model of the Edwards-Anderson model, on random graphs with arbitary 
degree distributions. The results are shown on the Nishimori line. We also show the results 
for the infinite-range model. 

§1. Introduction 

The study of spin models on complex networks has been carried out.'^^^ We 
study a spin model on random graphs with arbitary degree distributions as an ex- 
ample of the study of a spin model on complex networks. The behavior of spins on 
a no growing network is investigated. 

We investigate a Potts gauge glass modeP as a spin model. The Potts gauge 
glass model is a spin-glass model and is an extended model of the Edwards- Anderson 
model^' which is known as a spin-glass model. The understanding of the Edwards- 
Anderson model on random graphs and on the Bethe lattice is still incompleted .'^''^''^ 
The understanding of the Potts gauge glass model on random graphs and on the 
Bethe lattice is also incompleted. 

The Nishimori line-' is a line on the phase diagram for the exchange interactions 
and the temperature. The internal energy and the upper bound of the specific heat 
are exactly calculated on the Nishimori line.® The location of the multicritical point 
on the square lattice was conjectured, and it was shown that the conjectured value 
is in good agreement with the other numerical estimates. We will obtain results 
on the Nishimori line. 

There is a case where a percolation transition of networks occures. A network 
is divided into many networks by deleting some of its nodes and/or links. There is 
also a case where a percolation transition of clusters occurs. A cluster consists of 
fictitious bonds. The bond is put between spins. One of the clusters becomes a giant 
component when a cluster is percolated. We discuss the percolation transition of a 
cluster. 

We investigate the percolation transition of the Fortuin-Kasteleyn (FK) cluster 
in the FK random cluster modeL^'^^^ In the ferromagnetic spin model, the perco- 
lation transition point of the FK cluster agrees with the phase transition point. For 
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example, the agreement in the ferromagnetic Ising model is described in Ref. [TT]) . 
On the other hand, in the Edwards- Anderson model that has a conflict in the inter- 
actions, the percolation transition point of the FK cluster disagrees with the phase 
transition point. It was pointed out by de Arcangelis et al. that, despite the disagree- 
ment, the correct understanding of the percolation phenomenon of the FK cluster 
in the Edwards- Anderson model is important since a dynamical transition, which is 
characterized by a parameter called the Hamming distance or damage, is occurred at 
a temperature very close to the percolation temperature, and the dynamical transi- 
tion and the percolation transition are related to a transition for a signal propagating 
between spins.'^ We analytically obtain the percolation thresholds of the FK cluster 
for the Potts gauge glass model. 

We use a gauge transformation for deriving results. The gauge transformation 
was proposed in Ref. In addition to the application of the gauge transformation, 
results are shown by applying a criterioi]^ for spin models on the random graphs 
with arbitary degree distributions. 

In Ref. I13p . by applying the criterion with a gauge transformation, the percola- 
tion thresholds of the FK cluster for the Edwards- Anderson model on the random 
graphs with arbitary degree distributions were analytically calculated on the Nishi- 
mori line. 

We also show the results for the infinite-range model. 

This article is organized as follows. First in ^ a complex network model and 
the Potts gauge glass model are described. The FK cluster is described in ^ and 
appendix. A criterion for percolation of cluster is explained in 21 We will find in ^ 
the percolation thresholds. This article is summarized in ^ 

§2. A complex network model and a Potts gauge glass model 

A network consists of nodes and links. A link connected between nodes. The 
complex network model that we investigate is random graphs with arbitary degree 
distributions. The networks have no correlation between nodes. The node degree, k, 
is given with a distribution p{k). The links are randomly connected between nodes. 

We define a variable b{i,j), where b{i,j) gives one when nodes i and j is con- 
nected by a link. b{i,j) gives zero when nodes i and j are not connected by the link. 
The degree k{i) of node i is given by 



The coordination number (the average of the node degree for links), {k)^, is given 




(2-1) 



by 




(2-2) 



where Oat is the average over the entire network. N is the number of nodes. The 
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Fig. 1. Relation between the aspect a and the model on the network. 

average of the square of the node degree for Unks, {k'^)N, is given by 

1 ^ 

{k')N = ^Y.k\^). (2-3) 



N 

N 

We define 



2{k) 



N 



(2-4) 



{k^)N ' 

where a represents an aspect of the network. 

Figure [1] shows the relation between the aspect a and the model on the network. 
The network is almost a complete graph when a is close to zero, and the model on 
the network is almost an infinite-range model. The model on the network is the 
infinite-range model when {k)]\j = N — 1, {k'^)N = (-^ — 1)^, and a = 2/{N — 1). 
The network consists of many cycle graphs when a is one and {k)]\f is two. The 
model on the network consists of many chain models when a is one and {k) jsf is 
two. In the Erdos-Renyi (ER) random graph model and in the Gilbert model, the 
distribution of node degree is the Poisson distribution.'^' The ER random graph 
model is a network model wherein the network consists of a fixed number of nodes 
and a fixed number of links, and the links are randomly connected between the nodes. 
The Gilber model is a network model wherein the link between nodes is connected 
with a given probability. In the ER random graph model and in the Gilbert model, 
{k'')N = {k)N{{k)N + 1), and a = 2/{{k)N + 1). 

The Hamiltonian for a Potts gauge glass model, Ti, is given bj® 

^ = -- E E E e'f(-».^+«-'^>- , (2-5) 

i {i|b(i,i)=l}'^«.:;=l 

where qi denotes the state of the spin on node i, and qi = 0,1, q — 1. Uij denotes 
a variable related to the strength of the exchange interaction between the spins on 
nodes i and j, and Uij = 0, 1, . . . , g — 1. q is the total number of states that a spin 
takes. 

We use representations: Aj = e « and j'^ = Je i ^'^ Then, the Hamil- 
tonian (Eq. (j2-5p ) is given by 

, N q-l 

^ = -^E E E ^J'^^^r'^^P'^- (2-6) 

i {Mi,j)=i} n,j=i 
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The value of i^j.j is given with a distribution P{i'ij). The distribution P{uij) is 
given by 

P{^^,J) =pK,,o + - ^-..,o) , (2-7) 

where p is the probability that the exchange interaction between the spins is ferro- 
magnetic. 6 is the Kronecker delta. The normalization of P{vij) is given by 

q-l 

^P{uij) = l. (2-8) 

When Uij = {Jyj'" = J) for ah pairs, the model becomes the ferromag- 
netic Potts model. When q = 2, the model becomes the Edwards-Anderson model 
and is especially called the it J Ising model. 

In Ref. [T^ . it was pointed out that a gauge transformation has no effect on 
thermodynamic quantities. To calculating thermodynamic quantities, a gauge trans- 
formation wherein the transformation is performed by 

47^) ^ A,^Vi (2-9) 

— a- 

is used, where /ij = e <? , and qi is an arbitary valuve for qi. This gauge transfor- 
mation was proposed in Ref. Uj). By the gauge transformation, the Hamiltonian Ti 
part becomes V. ^ H. 

By using Eqs. ([23]), ([23]), and ([240]), the distribution P(fi,j) is given bjl^ 

P{i^ij) = Ae " ^'■.,.=1'''^ ^'^'■^^ , (2-10) 
where A and j3p are respectively 

^=^777—^ ^> (2-11) 



/3pJ, , ,8pJ ' 

e-?-l»-i) + (, - l)e— ^ 
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By performing the gauge transformation, the distribution P{vij) part becomes 
Yl P{T^i,j) = Ae " ^^^^=1 '-^^ 

{i,3) 

{Mi} 

where {x, y) denotes the nearest neighbor pairs conneced by links 



(2-12) 



JZ. g~^<i.J>^'-i,j=i-^i,i t^i ^ (2-13) 
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Fig. 2. Network and FK cluster. Three nodes, six links, three spins, an FK bond, and an FK 
cluster are depicted. Spins are aligned on each node and are represented by spin states. In this 
picture, the states of two spins are one and the state of a spin is zero. 

§3. The Fortuin-Kasteleyn cluster 

The bond for the FK cluster is put between spins with probabihty PpKiQi, Qj, '^i.j)- 
The value of PpK depends on the interaction between spins and the states of spins. 
We call the bond the FK bond in this article. PpKiQijQjii^ij) is given by 

ftK«„. = 1 - "•>^•"'•■''*«-^| , (34) 

where /3 is the inverse temperature, and /3 = 1/kBT. ks is the Boltzmann constant, 
and T is the temperature. By connecting the FK bonds, the FK clusters are gener- 
ated. In appendix, we will derive Eq. ()3-ip . By the gauge transformation, the Pfk 
part becomes PpK Pfk- 

Figure [2] shows a conceptual diagram of a network and an FK cluster. Three 
nodes, six links, three spins, an FK bond, and an FK cluster are depicted. Spins are 
aligned on each node and are represented by spin states. 

The thermodynamic quantity of the FK bond put between the spins on nodes i 
and j, [{bFK{iJ))T]R, is given by 

[{bFK{i,j))T]R = [{PFK{qi,qj,i^i,j))T]R, (3-2) 

where ( )t is the thermal average, and []r is the random configuration average. The 
thermodynamic quantity of the node degree for FK bonds at node i, [(/cFK(^))r]i?) 
is given by 

[{kFK{i))T]R = [{ ^ PFK{qi,qj,i'i,j))T]R- (3-3) 

The thermodynamic quantity of the square of the node degree for FK bonds at node 
h [(^fk(^))t]r, is given by 

[{k^FKii^hU 

J2 J2 PFK{quqj,J^i,j)PFK{qi,qi,i^i,i){'i- - Sj^i) 
{Mid)=i} mi,i)=i} 
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+ ^ PFK{qi,qj,i^i,j))T]R- (3-4) 
{i|fe(iJ)=i} 

The thermodynamic quantity of the node degree for FK bonds, [(A;fk)t]ro is given 
by 



1 ^ 

[{kFK)T]R = J^Y.^{kFKii))T]R. (3-5) 



The thermodynamic quantity of the square of the node degree for FK bonds, [(^fk)^]/?) 
is given by 

1 ^ 

[{4K)T]R = J^^[{4Kii))T]R. (3-6) 
i 

§4. A criterion for percolation of cluster 

We use a conjectured criterion for deriving the percolation thresholds. This 
criterion is a criterion of the percolation of cluster for spin models on the random 
graphs with arbitary degree distributions, and is given 

[{kLnd)T]R>2[{hond)T]R, (4-1) 

where /cbond is a quantity for a bond put between spins. A^bond for the FK bond 
is kpK for example. Equation ()4-ip is given by the inequality when the cluster 
is percolated. Equation ()4-ip is given by the equality when the cluster is at the 
percolation transition point. 

In addition, Eq. (|4Tp is true for sufficiently large number of nodes in the case 
that the magnitude of the bond does not depend on the degree k{i). 

We consider the condition that the magnitude of the bond does not depend on 
the degree k{i). We define a variable for the inverse temperature /3 as p{f3). We set 

0<p(/3)<l. (4-2) 

We consider a case that [(6bond(«, j))T]ij, [(A;bond(«))T]i?, and [(^bond(^))r]i? are re- 
spectively written as 

[{bbond{hj))T]R = p{^), (4-3) 

[{kbond{i))T]R = p{P)kir), (4-4) 

[{kLnAi))T]R = p\f3)k{imi)-l] 

+ p{P)k{i). (4-5) 

In the case, it is implied that the bias for k{i) does not appear in the statistical results 
of the bonds. Therefore, we describe the case that [(6bond(^) i))T]_R; [(^bond(^))T]R) 
and [(^bond(^))r]-R are respectively written as Eqs. ()4-3p . ()4-4p . and ()4-5p as the case 
that the magnitude of the bond does not depend on k{i). 

This criterion is a conjectured criterion and is not exactly derived yet. On the 
other hand, it was confirmecP^ that this criterion is exact for several extremal points 
when applied to the Edward-Anderson model. In this article, we do not examine 
this criterion and just apply this criterion to the present system. 
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§5. Results 

We will obtain the percolation thresholds of the FK cluster in this section. 

By using Eqs. ([231), (|2T3l) . and ([321), when /3 = (3p, the thermodynamic 

quantity of the FK bond put between the spins on nodes i and j, [(6fk(^> j))t]/?i is 
obtained as 



-/3pW({9l},{!^l,m}) 



{I'l.m} {Si} 

where Nb is the number of all links, and Nb = N{k) n /2. By using Eqs. dl^D, ([2T3l) . 
()3-ip . and ()3-3p . when /3 = /3p, the thermodynamic quantity of the node degree for 
FK bonds at node i, [(A;FK(^))r]/?5 is obtained as 

PFK«)Tk=(^^^^)A.W. (5-2) 

By using Eqs. (j2-9|) . (|2-13|) . (|3-1|) . and (|3-4|) . when (3 = /3p, the thermodynamic 
quantity of the square of the node degree for FK bonds at node i, [(^FK(^))r]-R; is 
obtained as 



e 



PpJ _l \ 2 



efpJ - 1 
e/3p-/ + q-l 



k{i)[k{i)-l] 

+ rU«. (5-3) 



When we set 



Eqs (|5-ip . (j5-2p . and (j5-3p are respectively formulated as Eqs. (j4-3p . ()4-4p . and 
(|4-5p . Therefore, the magnitude of the bond does not depend on the k{{). By 
using Eqs. ([331), dM]), (SiH), (E^D, and ([531), we obtain 



2(eiSpJ-l) + , - ■ 

Equation ()5-5p is given by the inequality when the cluster is percolated. Equa- 
tion (|5-5p is given by the equality when the cluster is at the percolation transition 
point. 
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From Eqs. ()4-2p and ()5-4p . there is the percolation transition point for < /3p < 
oo. From Eq. (j5-5p . there is the percolation transition point for < a < 1. By 
using Eqs. (|2-12p and (|5-5p . the probability p that the interaction is ferromagnetic is 
obtained as 

2(1 - a) + aq 

(2-a)q ^^^^ 
at the percolation transition point. By using Eqs. (|2-12p and (|5-6p . the percolation 
transition temperature Tp is obtained as 

Tp = ^ . (5-7) 



A;r In 



2(l-a)+ag 
2(l-a) 



Figure [3] shows the percolation thresholds of the FK cluster for the Potts gauge 
glass model. Figure [3l^a) shows the relation between the aspect a and the probability 
p. Eq. (|5-6p is used for showing Fig. O^a). Figure O^b) shows the relation between 
the aspect a and the percolation transition temperature Tp. Eq. ()5-7p is used for 
showing Fig. [3)^b). The solid line shows the result oi q = 3, the dotted line shows 
the result of g = 10, and the short-dashed line shows the result oi q = 100. J/kp is 
set to 1. 

For the ferromagnetic Potts model on the same network, the phase transition 
temperature 7"^^^™) i^E) 



rp(Ferro) 



kp In 



^ (5-8) 



2{l-a)+aq 
2{l-a) 



Tp (Eq. ()5-7p ) coincides with T^^^^°\ 

The complete graph is considerd as a ~ 0. We set {k)]\f = N — 1, {k'^)N 
{N — 1)'^, a = 2/{N — 1), and J —?- J/^/N. From the settings, the model on the 
network becomes the infinite-range model. By using Eq. ()5-6p . the probability 
that the interaction is ferromagnetic is obtained as 

,m . ^_il±i ^ 1 ,5.<,) 

(A^ - 2)q q 

for a sufficiently large number of nodes at the percolation transition point. By using 
Eq. (|5-7p . the percolation transition temperature Tp^^ is obtained as 

T^^^) = , (5.10) 

^ kpy/NHl + j^) kpq 

for a sufficiently large number of nodes. 

§6. Summary 



In this article, the Potts gauge glass model on the random graphs with artibary 
degree distributions was investigated. 
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The value of [(6FK(«,i)>r]R, [{kFK{i))T]R, [(^FK(^))r]i?, [(^FK)T]i?, and [(A;|,j^)r]i? 
on the Nishimori line were shown. They are quantities for the FK bonds and are 
exact even on a finite number of nodes. 

It is known that the internal energy and the upper bound of the specific heat 
are exactly calculated on the Nishimori line in the Potts gauge glass model without 
the dependence of the network (lattice).'^ In this article, it was realized that, as a 
property for the Nishimori line, the magnitude of the FK bond does not depend on 
the degree k{i). 

In this article, we showed the percolation thresholds of the FK cluster. It was 
shown that the percolation transition temperature Tp (Eq. ()5-7p ) on the Nishimori 
line for the Potts gauge glass model on the present network coincides with the phase 
transition temperature J'^^™™^!^) ^^j. ^.j^g ferromagnetic Potts model on the same 
network. The percolation thresholds of the FK cluster for the infinite-range model 
were also shown. 

We used a conjectured criterion Eq. (|4-ip to obtain the percolation thresholds. 
For this criterion, it was confirmed^"^' that this criterion is exact for several extremal 
points when applied to the Edward-Anderson model. Therefore, our entire set of 
results may be exact. 
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Appendix: the probabilities for connecting spins 



We will derive Eq. ()3-ip according to the method of Kawashima and Gubernatis 
in Ref. [T6|) . We define the weight of two spins on nodes connected by a link as 
w{qi,qj,iyij). w{qi,qj,iyij) is given by 



w{qi , qj , J- ) = exp I — 2^ 



exp 



27ri 



3-11 



We define the weight w for f jj- + qi — qj = as li^para- We obtain 

\ (3J{q-iy 
Wpi,ra[qi,qj,Uij) = ex.p . 

We define the weight w for + qi — qj ^ as i^anti- We obtain 

Wanti{qi,qj,T^i,j) = exp^-^^ . 



(6-3) 



We define the weight of graph for connecting two spins as ti'(5'conn)- We define the 
weight of graph for disconnecting two spins as w{gdisc)- We are able to write 



) + w{gdisc) , 
Wa.nti{qi, Qj^t^ij) = W(fi'disc) • 



(6-4) 
(6-5) 
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By using Eqs. ([621), (l631), dSlD, and ([621), we obtain 

'PJ{q-l) 



w{go 



exp 



exp 



q 



(6-6) 
(6-7) 



We define the probability of connecting two spins for i^ij + qi — qj = as i'paraC^conn)- 
We define the probability of connecting two spins for Vij + qi — qj 7^ as -Panti(5conn)- 
We are able to write 



-fpara (,9c 
-fanti (^c 



w{go 



w{g 

conn 

0. 



(6-8) 
(6-9) 



By using Eqs. (USD, dS^I), and ([63]), we derive Eq. (I3TD . 
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g. 3. Percolation thresholds of the FK cluster for the Potts gauge glass model, (a) The relation 
between the aspect a and the probability p is shown, (b) The relation between the aspect a and 
the percolation transition temperature Tp is shown. The solid line shows the result of g = 3, 
the dotted line shows the result oi q = 10, and the short-dashed line shows the result of g = 100. 
J/kB is set to 1. 



